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ABSTRACT 


One common feature of several pole placement techniques is the use of 
a dyadic (Rank 1) feedback matrix. The limitation of this design is exam- 
ined and a design involving output feedback matrices of Rank greater than 
on is developed as a logical extension of the dyadic feedback design. An 
f .ample is presented to illustrate the design procedure. 


INTRODUCTION 


Modal control has been suggested as a design tool in an effort to 
circumvent some of the problems in the design of multivariable systems using 
optimal control (Rosenbrock, 1962). In modal control the eigenvalues of the 
system matrix are changed to achieve the desii‘ed control objective. Much of 
the work which relates to the concept of modal control has been termed "pole 
placement." The design of linear multivariable control systems, with pole 
placement as a specification, has attracted the attention of several authors 
(Davison, 1970; Retallack and McFarlane, 1970; Fallside and Seraji, 1971; 
Sridhar and Lindorff, 1973). One common feature of these design methods is 
a dyadic feedback matrix. This paper discusses some of the disadvantages 
of this approach and extends the design procedure to feedback matrices of 
rank greater than one. 


STATEMENT OF THE PROBLEM 


Consider a controllable and observable linear time invariant multivariable 
system 


* •'*' A 4 

x = A x + Bu | 

(1 

y = Cx ) 

where x. is an n vector of states, u is an m vector f inputs and is 
a p vector of outputs. Let (X 1 ,A 2) . . -,A n ) and {pj,p 2 ,. • • ,Pn) the 
poles of the open loop and closed loop system, respectively. The problem of 
pole placement is to find the feedback matrix K such that the closed loop 
system matrix (A-BKC) has the eigenvalues (pi,P 2 ,. . . ,P n )- 
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CHARACTERISTIC POLYNOMIALS 


Let T be the n»n non-singular matrix of eigenvectors. We have 
open loop characteristic polynomial = |sl - A[ 


* (s - Xj)(s - X 2 ). . .(s - X n ) (2) 


and 


closed loop characteristic polynomial - | s I - A + BKCj 

= (s - pj) (s - p 2 ) . . - (s - p n ) (3) 


Further, 


[si - A 1 o [T* 1 1 |st-A| |T| = | s I - T~ 1 AT | = |sI-A| (4) 


where 


Also, 


A = diag 1X^X2,. . . ,X n ) 


| si - A+ BKC [ * | T“ 1 1 | si - A + BKC| ]T| = jf 1 (si - A + BKC)T| = | si - A + BKC| (S) 

where C = CT and B = T _1 B. Define M^BKC. Then, it is shown in the 
appendix that 
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where Mj , Mj j , Mjj^ . . . are the determinant of the matrices obtained by 
deleting the ith row and jth column, i,j rows and columns, i,j,k rows 
and columns, .... respectively. In addition, nu is the ijth element 
of M, J 


m i i 

m i j 

m i j 

m jj 


and A * (s - X j ) ( s - Xj 1* • *vS*^jj) 


For clarity, further development is carried out 'n terms of a system with 
distinct open loop eigenvalues. The open loop and closed loop character- 
istic polynomials are related by the equation 


si - A + BKC 
| si - a| 

The value of eq depends on the closed loop eigenvalues (pj,. . .,p n ). 
From equations (6) and (7), we have 
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Note that | Mj = 0 except for m = p = n. 


DYADIC K MATRIX 


Results for this cis*’ tin be obtained easily as a special case of 
equation (8). When K hrs rack one, that is. 


K = fd T 


(9) 


mi j = biKCj - b jfd T cj = d T Cjbjf 


(101 
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ith row of 


where f = r f 3 , f 2 * * * -.fm3 T > d T * [di,d 2 ,. . .,dn], bj is the 
B and cj is the ith column of C. With this choice of K 


| M| ,Mj ,Mj[ j ,Mj . . . ,M^ j - 0 


Now equation (81 reduces to 



Then, 


( 11 ) 


( 12 ) 


(13) 


Equation (13) is the same as equation (9) in Sridhar and Lindorff (1973) 
and K can be found following the procedure outlined by them. 

It is important to recall that the simplicity of the above pole- 
shifting algorithm is a consequence of K being a dyadic. However, this 
results in a loss of design freedom available in multivariable design prob- 
lems. In other words, K = fd T maps any output vector £ into a vector 
proportional to f . As £ varies the range of possible controls varies 
along the line f instead of the range of This is illustrated for 

m = p = 2 in figure 1. The effect of this is that the feedback will couple 
all the modes of the sys'/em although some of the modes of the system may be 
open loop decoupled and it may be difficult to satisfy performance criteria 
other than a simple allocation of the closed loop poles. In an effort to 
overcome some of these difficulties the design is extended to feedback 
matrices of higher rank in the next section. 


FEEDBACK MATRIX OF RANK 2 


In this case the determinants of all l * l matrices formed from M 
with l > 2 are equal to zero and equation (8) reduces to 


Next, 


± = M + f 




li 


(s - A j) (s - X j ) 


n m- 

+ s - A- 
i=l s A i 


M = 

A 


t 

1-1 


1 . J«l_ 

n s - Ai 

n (A a - a.) 
l=i J 


(14) 
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We have 


n-i n 
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( 16 ) 


By expanding the summation on the R.H.S. of (16) and regrouping the terms, 

H, 
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observing the fact that 


Mij = Mji 


mii m ij 


m j i mj j 


equation (17) becomes 
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Combining equations (14), (15), and (18) we have 
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Equating oefficients of (s - X-[) , 
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For n > 2, equation (19) reduces to 


Ax-Tx7 +m ii’ i = 1 ' 2 ’ 

1 = 1 1 J 


K can bn decomposed as the product of two matrices C and H, that is, 

K = GH where G is a (.m * 2) matrix and H is a (2 * p) matrix. Further, 


C = , H - [hj hj • • • hp] > JLi = [&il S i 2 1 » ll 


1 b 


6 



11*2 • * • JLlllp] 


illll 


K = 




&n!l2 * • • 


and m^j * bjKci = bjGHcj . Now, we have to solve the n nonlinear 
equations (20) in 2(m + p) variables for pole placement using output 
feedback matrices of Rank 2. 


Example 

This simple example illustrates the differences between the dyadic and 
Rank 2 method of pole placement. Consider the system 



The system transfer function is 


r i 


s 


1 


£ = 


0 


u 


1 2s ♦ 4 

s + 3 (s + 1) (s + 3) J 


The system is observable and controllable with m * 2, p = 2, and n = 3. The 
output y. depends only on the input The open loop poles are at 1, 

-1, and -3. It is desired to stabilize the system while maintaining that y^ 
depends only on uj . 

Max (m,p) « 2 poles can be placed using the output feedback. Let the 
closed loop poles be at -Pj = -0.5, -p 2 = -1.5 while the location of the 
third pole -p 3 depends on (-pj,-p 2 ). For this choice of closed loop poles, 

“l • M o * p 3 ) 
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For a dyadic feedback 


2 

u 2 * 32 + 03) 

°3 * 32 ( ' 3 P 3 ) 


Since 


We have 


. V 

f 2 J 


■ [$; i\t 


d ^c-jbj f = 04 , i = 1,2,3 


djfi =-3/2(1 + p 3 ) 
d 2 f 2 = 12/5(2 + p 3 l 
d 2 (f 1 + f 2 ) = l/10(p 3 - 3) 

Solving these equations with dj * 1 and d 2 = -2, we get 

f x = 0.6104 and f 2 = 0.0218 

For this choice of 

0.6104 -1.2208 

0.0218 -0.0436 



two poles can be placed at (-0.5, -1.5) and the third pole is located at 
-0.3023. However, it is not possible to satisfy the second criterion that 
yi depend only on ui as can bfe seen from the closed loop system matrix 


A - BKC = 


1 - fidj 
-f 2 di 

Cf 1 + 


-fid 2 
-1 - f 2 d 2 
- (fl ♦ f 2 )d 2 


-fjd 2 

-f 2 d 2 

-3- (f x + f 2 )d 2- 
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Next consider K with Rank 2 which can be written as 


pn f i 2 *|r a n ^12 ~| rin 

L f 21 ^22 J I ^2 1 ^22 J Lj&J 


[hi lL 2 l 


where 

ll 2 = [dj 2 
given by 


[fn f 12 )» jU = t f 21 ^22 1 * hi * [J\l d 2 i] T and 
T 

d 2 2 1 • The closed loop system matrix with the new feedback is 


A - BKC = 


1 - £ihi 

-izh. 


“£lll 2 

-1 * £ 2 ^ 




b(il + £2)hl -(&l + i2}h2 -3-(&i+a2%J 


If yi should depend only on Ui then 

-S.lll 2 “ 0 

From (20), for pole placement we have to satisfy the equations 
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These equations with the condition -g^h^ = 0 reduce to 

fllhl * &2II2 ♦ ilhl = “l 
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As in the dyadic case, the choice of closed loop poles at (-0.5, -1.5) 
require that 04 = (15/32) (1 + p 3 ), a 2 = (2/32) (2 + P->) and a 3 = (15/32) (-3 + p 3 ). 
Noticing the fact that £ihj = k n , j^h^ = k 22 anv - ki 2 , the feedback 
matrix has to be chosen subject to 
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K = 


and the third pole is located at -4. With k 2 i = -i/2, the closed loop 
system transfer function is 


X s 


1 

s + O.S 


s 2 + 1.55 - 5.5 


2s + 7 


(s + 0.51 (s + 1.5) (s + 4) (s + 1.5)(s + 41J 


u 


and the closed loop system meets the design requirements. 


OTHER METHODS 


Feedback matrices of higher rank can be easily computed for the pole 
placement of multivariable systems using state feedback by Kalman's (1971) 
method. The given system is transformed to the Luenberger canonical form 
(Jordan and Sridhar, 1973). Let Oj, o 2 , . . a m be the controllability 
indices of the system. Then the transformed system is given by 


l - A t + Bu 


where 
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The dimensions of Aj; and Bjj and (oj * oj ) and * m) , respectively. 
Further, 


and 



hi 


0 0 ... 0 
* 

0 ... 1 * * 

+ ith column 


The * denotes a nonzero element. From a study of the Luenberger's canonical 
form it is evident that the feedback affects only the oi,(oi+oV),. . ., 

(oi + 02 + . . o m ) rows of the system matrix. Let A m , B m , and A d consist 
of the oj,(oi + o 2 ),. . . ,(ox + 02 * • • • >o m ) rows of A, B, and (A + BK) , 
respectively. Then, 


A d ® A m + BjdK 


[ 21 ) 


and 
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where the lust row corresponds to the coefficient of the closed loop 
characteristic polynomial 


(s - P x ) (s - P2 1 . . . (S - t> n j * s n + a^s 11 * 1 + . . . UjS + a, 


and the l's appear in the 01 ,( 01 + 02 ) 
With this, A d , 


.,( 01 + 02 +. . . .ojji-j) columns. 


A + BK * 


0 1 . . . 0 

♦ » 
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o .... 1 
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has eigenvalues at ( 01 , 02 ,. • « >0n^ and K is given by (22] 
The above constructive procedure reduces to 


KC = B**(A d - A m ) 


(23) 


ter pole placement by output feedback. Munro (1973) has given conditions 
under which equation (23) can be solved to obtain arbitrary assignment of all 
the poles using constant output feedback. 

Sankaran (1974) has suggested an iterative procedure for determining the 
constant gain output feedback matrix that will stabilize a system. If the 
system matrix A varies to A + <$A, then the corresponding variation in the 
eigenvalue is given by (Rosenbrock, 1965) 


6 X j. 


trace[Q » 6A] 
trace[Q] 


(24) 


where 


n 

q = n 
i=l 
i#r 


(A - X.I) 


Assuming that <5A results due to output feedback and denoting the gain matrix 
by 6K, then 


6A * B • 6K • C 
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and a linear equation can be written for the elements of the gain matrix. 
Since the equation (24) is valid only for small perturbations, variations 
in the eigenvalues should be applied in small steps in the desired direction 
to compute 6K. 

The relative merits of the Jordan canonical form approach to pole 
placement outlined in this paper and the Kalman approach to pole placement 
are discussed in the next few paragraphs. Godbout (1974) has made an 
extensive comparison of the two methods. 

The Kalman method of pole placement uses the Luenberger Canonical Form 
of a system to shift the modes of that system. However, a system must be 
completely controllable in order to be transformed to Luenberger form. This 
restriction implies that the Kalman algorithm cannot be applied to any system 
that has at least one uncontrollable mode. This is not the case with the 
Joruan Canonical Form approach to pole placement. This algorithm allows the 
shifting of any or all of the controllable modes of a system that may also 
contain uncontrollable modes. This is a significant advantage for the 
Jordan Canonical Form approach over the Kalman approach. 

The control over pole motion is probably the most desirable feature of 
the Jordan Canonical Form scheme. The method allows the designer to choose 
the closed-loop value to which a specified open-loop mode will be shifted. 

The Kalman algorithm does not have this feature. The pole motion cannot be 
chosen or even determined when using this method. 

The Kalman pole placement scheme does have one desirable charac- eristic 
in that it is very straightforward. Furthermore, the Kalman method does not 
require Knowledge of the eigenvalues of the system, whereas the Jordan 
Canonical Form approach needs an accurate estimate of the eigenvalues in 
order to construct the Jordan form of the system properly. These eigenvalue 
estimates are generated by Francis' QR method which sometimes has difficulty 
converg nj, when a system has repeated roots. 


CONCLUSIONS 


This paper considers the problem of pole placement in multivariable 
systems using output feedback matrices of Rank greater than one with partic- 
ular reference to matrices of Rank 2. The results have been developed for a 
system with distinct eigenvalues. However, they can be easily extended to 
systems with multiple eigenvalues. The feedback matrices of higher rank 
provide considerable design freedom. But, the resulting equations for the 
feedback matrix are more complex than in the dyadic case. 
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APPENDIX 


Theorem: let ft he a diagonal matrix with elements (fi^.ftj,. . . ,ft n ) 

and let M be a n * n matrix. Then 
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where Mi, M^, Mjjfc . . . are the determinant of the matrices obtained by 
deleting the ith row and column, i,j rows and columns, i,j,k rows and 
columns, ...» respectively. In addition. 
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and A ~ ^1^2 • • * ^n* 

Proof: The proof depends on the repeated application of the following 

Lemma. 

Leima (Stiab, 1969) 

Let P, Q, and R be three n * n matrices identical except for their pth 
rows. Further, Lee the pth row of R be the vector sum of the pth rows 
of P and Q. Then, 
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that is, 

M 1 = |M| + ftjM! 

Again from equations (Al) and (A2) 
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Next , 
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Applying the lemma successively to M' 4 , M 5 , . . . , M n 
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Recalling the definition of ^ , 
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This proves the theorem. 

Closed Loop Chavact eristic Polynomial: 

| si - A + BKG J 

Th i s is a special case of the theorem with 

flj = (s - Xj), i = 1,2,. . .,n and M = BKC 
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